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Abstract 

We present a model describing a single brane with tension embedded into a five- 
dimensional space-time with compact extra dimension, which can be easily stabilized. 
We examine the linearized gravity in the model and obtain an expression for the four- 
dimensional Planck mass on the brane in terms of the model parameters. It is also shown 
that the scalar sector of the effective four-dimensional theory contain a tachyonic mode, 
and we discuss the problem of stability of the model. 

1 Introduction 

Brane world models and their phenomenology have been widely discussed in the last years. One 
of the most interesting brane world models is the Randall- Sun drum model with two branes, 
- the RSI model [1]. This model solves the hierarchy problem due to the warp factor in the 
metric and predicts an interesting new physics in the TeV range of energies. 

Most of the brane world models with one compact extra dimension and thin branes with 
tension demand the existence of at least two branes. At the same time the matter located 
on the brane, which is not "our" brane, can strongly affect the world located on "our" brane. 
For the case of the RSI model it was shown in [2]. So it would be quite interesting to find 
out, whether it is possible to construct a model with only one tensionful brane in a compact 
extra dimension, admitting a solution to the hierarchy problem in the way analogous to that 
proposed in [T]. 

A characteristic feature of models with single brane is the presence of at least one tachyonic 
mode in the perturbative linearized theory [3|. At the same time the linearized theory, as well 
as the five- dimensional effective action describing a brane world model, is valid for the energy 
range of the order of the fundamental energy scale of the theory, defined by the five- dimensional 
gravity (we suppose that this scale is of the order of 1 — 10 TeV). Thus, if the masses of the 
tachyonic modes are far beyond the energy range of its applicability, their influence on the 
theory cannot be accessed in the linear approximation, and one needs to consider the nonlinear 
effects. 

Some solutions with single brane in a compact extra dimension, interesting from the cosmo- 
logical point of view, were obtained in |3] . But the energy- momentum tensors used for obtaining 
these solution are " phenomeno logical" , i.e. they are added to the action "by hand". 

Here we present a model describing the scalar field minimally coupled to gravity in a five- 
dimensional space-time, admitting the existence of a single brane quite naturally and being 
of interest from the point of view of the hierarchy problem. Moreover, the size of the extra 
dimension in this model can be easily stabilized. Thus, the model appears to be devoid of 
the main fiaw of the original Randall-Sundrum model - the existence of the massless scalar 
mode, called the radion, which arises due to the fiuctuations of the branes with respect to each 
other and whose interactions contradict the existing experimental data. We argue that the 



stabilization of the size of extra dimension is made in the same way as in [5] . This method is 
free from the main disadvantage of the approach proposed in [H], where the backreaction of the 
scalar field on the background metric is not taken into account. There is only one tachyonic 
mode in the model with the mass of the order of the four-dimensional Planck mass ~ lO^^GeV, 
thus lying far beyond the applicability range of the theory. 

The paper is organized as follows. In Section 2 we present the background solution and the 
method of its stabilization. In Section 3 we obtain gauge conditions and equations of motion 
for linearized gravity in the model. In Section 4 we consider the tensor modes and obtain an 
expression for the four- dimensional Planck mass on the brane in terms of the model parameters. 
In Section 5 we consider the scalar sector of the theory, discuss the stability of the model and 
obtain the estimates for the mass of the lowest scalar mode and its coupling constant to matter 
on the brane. And finally, we discuss the obtained results. 



2 The model 

Let us denote the coordinates in five- dimensional space- time E = M^x S^hy {x^} = {x^',y}, 
N = 0,1,2,3,4, /i = 0,1,2,3, the coordinate = y,--L < y < L parametrizing the fifth 
dimension with identified points —L and L. The brane is located at the point y = L. 
The action of stabilized brane world model can be written as 



S = I d X I dy^/^ 



^ 1 

2M^R--g''''dM<pdM<p-V{<P) 



y=L 



Here V{(j)) is a bulk scalar field potential and A(0) is the brane scalar field potential, g = detg^y, 
and denotes the metric induced on the brane. The signature of the metric guN is chosen 
to be (-,+,+,+,+). 

The standard ansatz for the metric and the scalar field, which preserves the Poincare in- 
variance in any four- dimensional subspace y = const, looks like 

ds'^ = e-^^^y'>7]^ydx^'dx'' + dy^ = -fMN{y)dx^^dx^ , (2) 
(f){x,y) = (f){y), (3) 

7]^^ denoting the flat Minkowski metric. If one substitutes this ansatz into the equations 
corresponding to action ([T]), one gets a rather complicated system of nonlinear differential 
equations for functions A{y), (f){y): 

dV dX ^. ^, ^ ^, ,, ,,, 

i2M'{Ar+^iv-^{<pr) = o, 



\ (^(0')' + V + X6{y - L)^ = -2M' {-3A" + 6(A')^) . 



(4) 



Here and below ' = d/dy. An interesting conclusion following from these equations is that the 
relation 

holds in the bulk for any potential V{(j)), and thus A" > in the bulk. This inequality was also 
obtained in [7] from the weaker energy condition. 
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To find an analytic solution to this system we will use the results of O [8] . Let us consider 
a special class of potentials, which can be represented as 



Let us also suppose that 



w{<p) = (6) 



In this case the scalar field potential takes a simple polynomial form 



and the corresponding continuous background solution can be easily found (with the help of 
the procedure described in [5|,[H]): 

4>^byf-. (8) 

The additive constant in the solution for A{y) is chosen in such a way that the coordinates 
{x^} are Galilean on the brane (see [HI [ID] for details). We will refer all the energy parameters, 
which appear in the theory, to this Galilean coordinate system on the brane. 

In order the equations of motion be valid on the brane too, one needs to finetune the brane 
potential A ((/)). We choose 

A(0) = -Wi<P) = (9) 

In this case the brane appears to be of the BPS type. The size of the extra dimension is not 
defined by the solution yet. 

To stabilize the size of the extra dimension, let us add the following term to the scalar field 
potential on the brane: 

AA(0)=/32(0-0o)'. (10) 
Such an addition will not affect the equations of motion provided 

0|y=L = 00, (11) 

which means that 



L = ^. (12) 

7 

Thus, we see that the size of the extra dimension is stabilized. 

It is necessary to note that the background solution presented above was obtained without 
imposing Z2 orbifold symmetry, which is inherent to the most brane world models, although 
the solution itself possesses reflection symmetry with respect to the point y = 0. 

We also suppose that the parameters of the potentials 7, 0o, /3, when made dimensionless by 
the fundamental five- dimensional energy scale of the theory M, should be positive quantities 
of the order 0(1), i.e. there should be no hierarchical difference in the parameters. We note 
that action ([T]) and the corresponding four-dimensional effective theory can be used only at the 
energy scales E < M measured in Galilean coordinates on the brane. 
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3 Linearized gravity 



Now let us turn to the examination of linearized gravity in the model. We represent the metric 
and the scalar field as 

9MN{x,y) = -fMN{y) + ^y=ihMN{x,y), (13) 
<j){x,y) = <p{y) + ^=f{x,y). (14) 

To simplify the analysis, let us impose Z2 orbifold symmetry conditions (although this symmetry 
is not necessary for obtaining the background solution). Correspondingly, the metric qmn and 
the scalar field (f) satisfy the orbifold symmetry conditions 

g,iu{.x, -y) = g^^{x, y), g^^ix, -y) = -g^^ii^x, y), 

gu{x, -y) = guix, y), 0(x, -y) = 0(x, y). (15) 

We realize that imposing Z2 orbifold symmetry is a rather artificial procedure. But a 
consistent and thorough analysis of linearized gravity without this symmetry, i.e. taking into 
all the degrees of freedom coming from the metric, is a very complicated problem (for example, 
we cannot impose the gauge conditions which will be used later). At the same time, a theory 
with the orbifold symmetry makes sense and was studied, for example, in [3]. Moreover, we 
have a developed formalism for studying linearized gravity in brane world models stabilized by 
the bulk scalar field and with extra dimension forming the orbifold /Z2 - see [TT]. The only 
difference from this case is that all the fields should have a "good" behavior at the point y = 0, 
i.e. the fields should be smooth at ?/ = 0, which corresponds to the absence of the brane as a 
physical object at this point. For these reasons, in this paper we restrict ourselves to the case 
with Z2 orbifold symmetry conditions. 

Substituting representation ( fT3|) and into action ([T]) and keeping the terms of the second 
order inhuN and /, we get the second variation Lagrangian of this action pjj. This Lagrangian 
is invariant under the gauge transformations 

hMNi^^y) = hMN{x,y) - {Vm^N + Vat^m), 
f('\x,y) = f{x,y)-<P'U, 

where Vm is the covariant derivative with respect to the background metric, provided C,M{x,y) 
satisfy the orbifold symmetry conditions 

^ti{x,-y) = ^f,{x,y), U{x,-y) = -U{x,y). 

These gauge transformations are a generalization of the gauge transformations in the unsta- 
bilized RSI model P, CD]- We will use them to isolate the physical degrees of freedom of the 
fields Hmn and /. We also note that since ^4\y=L = 0, the brane appears to be straight (the 



disadvantages of bent-brane formalism were discussed in [12j). 

It was shown in [11] that with the help of these gauge transformations one can impose the 
gauge 

ie-'^h,,y - ^s^-'Wf = 0, (16) 

= 0, 
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after which there remain the gauge transformations satisfying 



(e^,)' = 0. (17) 

A substitution 

h^u = V ~ ^ItiuhiA (18) 

allows us to decouple the equations of motion, following from the second variation Lagrangian, in 
gauge ( fT6l) . Gauge transformations satisfying (ITTj) allow one to impose the traceless-transverse 
gauge condition on the field h^y [101 [H] 

h = l^'b,, = 0, d-'b,, = 0, (19) 

the residual gauge transformations now being 

= e~'^e^(a;), 9%(x) = 0, De. = 0, (20) 

where □ = 'q'^^d^du- Transformations ( l20i) act only on the massless mode of the field 6^,^ and 
provide the correct number of degrees of freedom of the massless graviton |lUj . 

Finally, we get the equations of motion in the interval (0, L) with corresponding boundary 
conditions at the points y = 0, y = L for the field b^^ 



(21) 



&;j.=+o = 0, (22) 
b'+2A\X=L-o = 0, 



and for the field g = e ^^*^^)/i44(x, ?/) 

6M3 



9" + 2g'{A'-^)-^^g + d,d^g = 0, (23) 



9'\y=+0 



/3V - d^d^g\y=L-o = 0, (24) 

see [TT] for details. 

4 Tensor modes and the hierarchy problem 

Let us study first the modes of the tensor field bfj_y{x,y), which satisfies Eq. (1^ . Substituting 
into this equation 



bf,u{x, y) = c^ue'^^'ipniy), Cf.^ = const, = -m. 



2 

m 
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we get: 

<|,=+o = < + 2AVn|,=L-0 = 0. (25) 

The boundary conditions suggest a substitution ipn = exp(— 2A)c<j„ (note that A'\y=+o = 0), 
which turns this equation into 

d 

dy 

uj'Jy=+o = ^n\y=L-o = 0. (26) 

We see that the eigenfunctions Un are solutions of a Sturm-Liouville problem with von Neumann 
boundary conditions. In accordance with the general theory [13j, the problem at hand has no 
negative eigenvalues for arbitrary A, only one zero eigenvalue, corresponding to uo = const. 

The eigenfunctions {"ipniu)} of eigenvalue problem build a complete orthonormal set, 
the eigenfunction of the zero mode being 

My) = Ne-'^^y\ (27) 

Expanding in this system 

oo 

b,. = Y,b;,{x)My), (28) 

n=0 

we get four- dimensional tensor fields b'^j^{x) with definite masses. 

A standard technique gives us an expression for the four- dimensional Planck mass on the 
brane 

rL f°° (7!/)'* 

Mpi = M^ e'^^dy ~ M^2e^ / e~T^dy= (29) 



J-L Jo 



47 v4/ 7 



and 



,7 



Mp; f« 2M e^^. (30) 

Let us suppose that all fundamental parameters of the theory lie in the TeV range. To have 
the hierarchy problem solved, one should take 

' - ^° -36, (31) 



36M3 36M3 
which means that 



36M2 (32) 
and 



~ . (33) 
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Although Eq. ( l26i) cannot be solved analytically for n 7^ 0, it is reasonable to suppose that the 
lowest masses of the four-dimensional tensor excitations b^j^{x) are of the order of L~^. 
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5 Scalar sector and stability 



In order to find the mass spectrum of the scalar particles described by Eq. ( |23l) let us substitute 

g{x,y) =e'P''gn{y), = -f^l, 

into this equation: 



.'\2 



(34) 



h^n^ yn\y 



9n\y=+0 '■ 
L-0 



0. 



(35) 
(36) 

0, from ([ISD it 



It is necessary to note that since the field / should be smooth at the point y 
follows that the value {g'^J (j)') should be continuous at y = too. 

First, let us solve Eq. flH^ for the case /xq = 0, i.e. for the zero mode. In the case of 
background solution ^ the wave function g^, satisfying boundary condition at ?/ = 0, has the 
form 

4 4 

■} „ 'I y 

"'dq. (37) 



3 \y\ 

14 



<1 
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(18M3)' 

It is not difficult to check that 5'o|y=L 7^ 0. Thus, the scalar zero mode is absent in the model. 

Now let us examine, whether there are scalar tachyons in the model. To this end we denote 
fi^ = — /i^ > (here and below we omit the subscript n) and introduce a new dimensionless 
variable 

t = — o-y. 

Mi 



(38) 
(39) 
(40) 



In this case Eq. (1341) and boundary conditions take the form 



9 + ^9 
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"^fg-fi^exp ]g = 0, 



-2 
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g + fi exp ( — 1^ 



0, 



7L 



Ml 

9\t=o = 0, 



where /i = /i^^ exp (^—^^^ , ~ 6 (see (!33|) ) and g = dg/dt. 

Unfortunately we cannot solve Eq. ( l38l) analytically. Numerical analysis (see Appendix A) 
shows that for Ji < 0.9507, t = 3.2: g{t) > and g{t) > (see examples on Figs. [H At the 
same time for fl > 0.9509, t = 3.2: g{t) < and g{t) < (see examples on Figs. El HI Fig. |l]is 
shown for t < 2, but one can check that for /2 = 1.5 and t = 3.2 g{t) < and g{t) < 0). The 
graphs on Figs. [H Ej [3] are shown for t < 3.3, it is made to show the behavior of g{t) in the 
interval t G [0,3.2], especially for the cases p, = 0.9507 and /2 = 0.9509. For t > 3.2: g/g > 0, 
it can be easily seen from the structure of Eq. fl38|) . Indeed, let us divide fl38|) by g and pass to 
the equation for q{t) = g/g, which takes the form of a Riccati equation 



q + q^ + 2 



-t^ + yU^ exp 



(41) 
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If initially for some value of t = tg > 0: q > 0, then q will remain positive for any t > tg. Indeed, 
the function q should pass through zero to change the sign. But if < g -C 1, then from Eq. 
f l4ip it follows that q > 0, q appears to be a growing function and thus remains positive. Thus, 
for t > 3.2: g{t) > 0, g{t) > or g{t) < 0, g{t) < depending on the sign of g{t) at t = 3.2. In 
both cases ( l39l) is not satisfied, since M > 0, 7 > and > 0. 

Eq. (!39l) can be satisfied if g{t) = g{t) = at the point t = ^ 6. But such boundary 

conditions imply that g{t) = in the interval t G e, for any e > - this conclusion follows 

from the theorem about the existence and uniqueness of solution for the Cauchy problem - see, 
for example, [U]. Finally, due to the continuity of the function g{t), we get g{t) = in the 

interval t G 0, ^ . Of course, the same conclusion can be made for any < ti < such 

that g{h) = g{h) = 0. 

As for the region 0.9507 < Ji < 0.9509, we have made a large number of numerical simula- 
tions with different values of Ji. The behavior of the corresponding solutions is such that for 
t > 3.6: g{t) > 0, g{t) > or g{t) < 0, g{t) < respectively, analogous to the behavior of 
solutions presented on Figs. [21 [3l A simple qualitative explanation of this fact can be given. 
For 0.9507 < Ji < 0.9509 and t > 3.6 the coefficient /i^exp {jg^ in (1381) grows rapidly, which 

leads to the growth of the absolute value of function g{t) with coordinate t for t > 3.6. 

Nevertheless, for some value of Ji such that 0.9507 < ft < 0.9509 there exists a solu- 

tion, which satisfies condition ( 1391) . Indeed, let us define a function F{p,) = g{t) + 

/i^exp [^]g{t)\i^ 7i . For ft < 0.9507 it is positive, whereas for 0.9509 < /i it is negative 

(see Figs. HI 12, 131 H]) . Thus, it is reasonable to suppose that there exists an appropriate value 
fi* (0.9507 < Ji* < 0.9509) such that F{fi*)\^^ ~,l = 0, which corresponds to a tachyonic mode. 

A/ 3 

It appears to be very difficult to find the exact value of the tachyonic mass numerically. At the 
same time the physical mass of the tachyon is such that 

/x*^ ^ ^ -^M^, ^ - (0.5 ■ lO^'GeV)' (42) 

7 

for the given values of the model parameters (we suppose that 7 ~ M*). Such energy scale 
lies outside the range of validity of our effective theory, described by action ([T]) (because ~ 
E ^ M, see Section 2). From the classical point of view it can be understood as follows: the 
tachyonic mode should behave as e'^*^ . The time derivative of the tachyon field ~ /i*e'^*^ , i.e. 
it is enhanced by the large value of /i* ~ Mpi in comparison with the tachyon field itself. Thus, 
even if the value of the tachyon field is small, its time derivative would lead to breakdown of 
perturbative approach and corresponding nonlinear effects, coming from the five- dimensional 
curvature (through substitution ( |T8l) ). Another remarkable thing is that the wave function of 

the tachyon is such that if g\y=o = 1, then in the leading order g\y=L ~ exp (^~/^* exp ^Ig^j j ~ 

exp (— /i* exp (36)). It means that when the nonlinear effects and (or) effects of the underlaying 
fundamental theory begin to affect the behavior of the theory in the bulk, the theory on the 
brane remain intact, because the coupling constant of the tachyon to matter on the brane, 
which is proportional to the value of the wave function on the brane, is negligibly small - much 
smaller than the coupling constant of the massless tensor graviton. Thus, the runaway of the 
scalar field can be stopped in the bulk because of the nonlinear effects coming from action ([T]) 
or from the underlaying fundamental theory. Of course, we cannot argue that it is indeed so. 
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0.5 1 1.5 2 2.5 

Figure 1: Numerical solution ior g{t), jl — 0.5 
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but such situation can be realized. 

Of course, our examination is not explicit since it is based on the numerical calculations. 
But we think that the analysis made testifies in favor of absence of scalar tachyons in the model 
below the energy scale of our effective theory ([T]). As for the ghosts, the form of the effective 
action for the scalar modes ensures the proper signs of the appropriate kinetic terms pjj. 

The form of Eq. (138!) allows us to estimate the mass of the lowest scalar excitation and its 
coupling to matter on the brane. Indeed, let us suppose that the lowest mass fii (see (Ell)) is 
such that Hi/M ^ 0{1). In this case we can neglect the last term in Eq. (I34p in comparison 
with the last but one term of this equation, and the solution of the resulting equation takes the 
form 



g^iy) ^A,\e ism'^ + / q te "rfg | ~ goiy), (43) 

(18M3)3 Jo 

where Ai is a normalization constant. Let us suppose that the size of the extra dimension is 
such that 7L/M4 = 6, see (!33l) . The values gi{L) and g[{L) can be easily calculated, which 
gives us 

gi{L) ^ A, ■ 89.Q, (44) 
g[iL)^A,-AA.8^. (45) 

Substituting (Hll) and (H5l) into (j36l) we easily get 



.2 



H-i — 3"- (4o) 

For example, if ~ M, 7 ~ M'^/^ and M ^ lOTeV, the lowest mass fxi ~ 7TeV. Of course, 
it can be even smaller depending on the values of the parameters (3, 7 and M. 

It is also necessary to note that the analysis carried out with the help of the numerical 
solution of Eq. flMl) for such small /ii reproduces the results obtained using fj^3|) with a very 
good accuracy (of the order of 1 — 2%). 

Now let us calculate the coupling constant of the first scalar mode to matter on the brane. 
To this end we need to calculate the normalization constant Ai. The normalization condition 
for the scalar modes takes the form fTT| 





It is more convenient to pass to the variable t = —^y: 

6 



[9i + 7JiS^9i = (47) 



3 





The integral in (HHll can be evaluated numerically, which gives us 

A? ^ 0.004^. (49) 

Now we can calculate the coupling constant of the lightest scalar mode to matter on the brane 
(see HH): 

„ = .J^ . . _,pl. (50) 
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One can see that ei ^ joYeV given values of the fundamental parameters 7 and M. 

Unfortunately it is impossible to calculate even the lowest mass of the tensor excitations 
using the method described above. One should carry out a very precise numerical analysis to 
get an information about the spectrum of the tensor modes. 

6 Conclusion 

In this paper a model describing the scalar field minimally coupled to gravity in the spacetime 
with one compact extra dimension is proposed. It admits the existence of a single tensionful 
brane, contrary to the most brane world models with one compact extra dimension demanding 
the existence of at least two branes (of course, except the simplest case of the ADD model 
P!H| [TU] with tensionless branes). We also showed that the model could be interesting in view 
of the hierarchy problem. 

The linearized gravity in the model was studied under the assumption of the Z2 orbifold 
symmetry. We obtained the expression for the four-dimensional Planck mass on the brane in 
terms of the fundamental five- dimensional parameters of the theory. We also made a stability 
analysis of the model, - analytical for the tensor modes and numerical for the scalar modes, 
which resulted in the conclusion that the scalar sector of the model contain one tachyon, which 
corresponds to the result obtained in [3], its "mass" being of the order of the four- dimensional 
Planck mass. Thus, the model as it is, at least in the linear approximation, is unstable and its 
" lifetime" is of the order of the four- dimensional Planck time. Nevertheless, the energy scale of 
the tachyon is such that multidimensional nonlinear fundamental underlaying theory can come 
to play and " lift up" the scalar sector from falling down. Of course it is not necessarily so, but 
in principle it seems to be possible. 

The background solution can also be used to describe the world with two branes. Indeed, 
the second brane can be placed at the point y = Lq, < Lq < L. The results of [H] suggest 
that such a system is totally devoid of tachyons. 

It is very interesting to carry out a numerical calculation of the coupling constants and 
the masses of the tensor modes and a complete description of the scalar sector of the model, 
as well as the model without Z2 orbifold symmetry (in the latter case there should appear 
antisymmetric modes). One can also use the model discussed in this paper (for example, the 
stable configuration with two branes) as a basis for constructing models with universal extra 
dimensions. These tasks deserve additional thorough investigation. 
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Appendix: note on numerical analysis 

We solve Eq. (!38l) numerically with the following initial conditions on the "time" variable t: 

git)\t=o = i, (51) 

m\t=o = 0. (52) 

But since the coefficient 

^ _1 
9 t 

in flHSl) is not defined at the point t = 0, it is inconvenient to use the point t = as the initial 
point for numerical calculations. To bypass this problem, we find an approximate analytical 
solution of Eq. (l38l) in the vicinity of the point t = 0: 

g{t) ^ 1 - P^t'. (53) 

Now we choose the point to = 10~^^ as the initial point instead of t = 0. The corresponding 
initial conditions take the form 

^(^)k=.o = l-ylO^'^ (54) 

g{t)\t=to = -^'lO-''- (55) 

The numerical analysis for a large number of different values of fl was made using the program 
package Mathematical version 5.2. Selected solutions are presented on Figs. [Il[2l[3llll 
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